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ABSTRACT 

We study in a systematic way the quality factor of the lower and upper kHz QPOs 
in a sample of low luminosity neutron star X-ray binaries, showing both QPOs vary- 
ing over a wide frequency range. The sample includes 4U 1636-536, 4U 1608-522, 
4U 1735-44, 4U 1728-34, 4U 1820-303 and 4U 0614+09. We find that all sources 
except 4U 0614+091 show evidence of a drop in the quality factor of their lower kHz 
QPOs at high frequency. For 4U 0614+091 only the rising part of the quality factor 
versus frequency curve has been sampled so far. At the same time, in all sources but 
4U 1728-34, the quality factor of the upper kilo- Hz QPO increases all the way to the 
highest detectable frequencies. We show that the high-frequency behaviours of both 
the lower and upper kHz QPO quality factors are consistent with what is expected if 
the drop is produced by the approach of an active oscillating region to the innermost 
stable circular orbit: the existence of which is a key feature of General Relativity in the 
strong field regime. Within this interpretation, our results imply gravitational masses 
around 2Mq for the neutron stars in those systems. 
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1 INTRODUCTION 

Kilohertz Quasi-Periodic Oscillations (kHz QPOs) have now 
been detected by the Rossi X-ray Timing Explorer (RXTE, 
Bradt et al. 1993) from about 25 neutron star low-mass 
X-ray binaries (van der Klis 2006 and references therein). 
These signals, whose timescales correspond to the dynam- 
ical times of the innermost regions of the accretion flow, 
have triggered much interest because they may carry im- 
prints of strong field general relativity, such as the existence 
of an innermost stable circular orbit (ISCO) around a suf- 
ficiently compact neutron star. For instance, a drop in the 
amplitude and quality factor of the QPOs at some limiting 
frequency was proposed as a possible signature of the ISCO 
(e.g., Miller, Lamb, & Psaltis 1998). 

Following this idea, we studied in a systematic way the 
QPO properties of 4U 1636-536, with emphasis on the de- 
pendency of the quality factor (Q = I//FWHM) and ampli- 
tude with frequency (Barret, Olive, & Miller 2005a,b). We 
have thus shown that the lower and upper kHz QPOs of 
4U 1636-536 follow two distinct tracks in a Q versus fre- 
quency plot. The quality factor of the lower kHz QPO in- 
creases with frequency up to a maximum of Q ~ 200 at 
~ 850 Hz, then drops precipitously to Q ~ 50 at the highest 
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detected frequencies Slower ~ 920 Hz (a similar behaviour 
was seen earlier from 4U1608-52, Barret et al. 2005c) On 
the other hand, the quality factor of the upper kHz QPO 
increases steadily all the way to the highest detectable QPO 
frequency, although the quality factor is lower than for the 
lower QPO. The rms amplitudes of both the upper and lower 
kHz QPOs decrease steadily towards higher frequencies. 

In this paper, we extend the analysis carried out for 
4U 1636-536 to all low-luminosity neutron star binaries for 
which kHz QPOs have been detected over a large range of 
frequencies. The sample includes: 4U 1636-536, 4U 1608- 
522, 4U 1820-303, 4U 1735-44, 4U 1728-34, 4U 0614+09. 
An exhaustive list of references on QPOs from these sources 
is available in van der Klis (2006). In § 2 we describe our 
analysis procedure, and in § 3 we discuss the implications of 
our results and show how the high-frequency dropoff in Q 
can be accommodated quantitatively in a toy model based 
on the approach to the ISCO. 



2 DATA ANALYSIS 

For the purpose of this paper, we have retrieved all science 
event files from the RXTE archives up to the end of 2004 for 
all six sources. The files are identified with their observation 
identifier (obs ID) following the RXTE convention. An Obs 
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Figure 1. Quality factor versus frequency of all fitted QPOs: red circles identify fitted lower QPOs, blue circles identify fitted upper 
QPOs, black circles correspond to single QPOs. Each data point represents the average over a continuous data segment. 



ID identifies a temporally contiguous collection of data from 
a single pointing. Only files with time resolution better than 
or equal to 250 micro-seconds and exposure times larger 
than 600 seconds are considered. No filtering on the raw 
data is performed, which means that all photons are used in 
the analysis. Type I X-ray bursts and data gaps are removed 
from the files. 

For each file identified with its Obs ID, we have 
computed Leahy normalized Fourier power density spectra 
(PDS) between 1 and 2048 Hz over 8 s intervals (with a 1 
Hz resolution). N 8-second PDS are thus computed. N is 
typically around 400 in most files, whose duration ~ 3200 is 
consistent with the orbital period of the RXTE spacecraft. 
A Fourier power spectrum averaging the N PDS is first com- 
puted. This averaged PDS is then searched for a QPO using 
a scanning technique which looks for peak excesses above 
the Poisson counting noise level (Boirin et al. 2000). No fit 
is performed at this stage as the scanning procedure returns 
the centroid frequency of each excess and an approxima- 
tion of its spread in frequency (a rough measure of its total 
width) . In case of the presence of two excesses, the one with 
the higher significance is considered: its centroid frequency 
is vo and its spread is wo- 

We now wish to attribute to all TV PDS the best esti- 



mate of the instantaneous QPO frequency v qp0l N to apply 
later on the shift-and-add technique. To start, u qpoi N is 
set to v>o- We define a window of 25 Hz on both sides of 
the excess centered at vq and recursively search for excesses 
over time intervals of shorter and shorter durations, i.e. in a 
smaller and smaller number of 8 second PDS averages. More 
specifically, in the first iteration, we thus consider two con- 
secutive time intervals averaging the 8-sec PDS from number 
1 to N/2 to produce PDSi,i and the 8-sec PDS from num- 
ber N/2 + 1 to N to produce PDSi, 2 . In PDS M and PDSi, 2 , 
we search for peak excesses between vq — wo/2 ~ 25 Hz and 
zA) + ™o/2 + 25 Hz, still above the Poisson counting noise level 
at the 4<t level. For instance, if an excess is detected PDSi,i 
at a frequency 1/1,1 (spread wi,i), then we set v qp0l N/2 to 
1/1,1 before the next iteration. If on the other hand, no excess 
is detected in PDSi,2, then v qp0¥j/2+l N remains unchanged 
at vq. In the next iteration, we repeat the procedure con- 
sidering two further intervals averaging the 8-sec PDS from 
number 1 to N/A to produce PDS2,i and the 8-sec PDS 
from number 7V/4+ 1 to N/2 to produce PDS2.2, and search 
in both PDS an excess between 1^1,1 — uii,i/2 — 25 Hz and 
vi t i + wi,i/2 + 25 Hz. The tracking procedure stops when 
no more excess is detected in any of the intervals analyzed. 
At the end, the procedure thus returns the best possible 
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estimate of the instantaneous QPO frequency in each 8 sec- 
ond PDS. We have checked through simulations of a QPO 
signal of varying frequency and of amplitude similar to the 
real data that this procedure follows with great accuracy the 
changes in frequency. 

In order to reduce the contribution of the long term fre- 
quency drift to the broadening of the QPO profile, within 
a data file, we then shift-and-add (Mendez et al. 1998) the 
individual 8 second PDS to the mean QPO frequency over 
the file. The shifted and added PDS is then searched for 
QPOs with the scanning routine. A quality check is per- 
formed at this stage with the scanning routine to verify that 
the tracked QPO in the resulting PDS has a centroid fre- 
quency consistent with uo and a spread smaller than wo- 
Either one or two QPOs are then fitted with a Lorentzian 
of three parameters (frequency, full width at half maximum, 
and normalization) to which a constant is added to account 
for the counting noise level (close to 2.0 in a Leahy nor- 
malized PDS). The fitting is carried with the XSPEC 11.3.2 
spectral package (Arnaud 1996), taking advantage of the ro- 
bustness of its fitting procedures (including the error com- 
putations). Next, we keep QPOs which are detected above 
2.5a (the significance being defined as the integral of the 
Lorentzian divided by its error), above 500 Hz, with a qual- 
ity factor larger than 3. 

In Figure 1, we show the quality factor of the QPOs 
detected for all six sources. Although with fewer details than 
for 4U 1636-536, the same trends are present for all the other 
sources as well. In particular a sharp drop of the quality 
factor of the lower kHz QPO is also seen in 4U 1608-522 and 
4U 1728-34. The drop is also seen in 4U 1735-44 and to a 
lesser extent in 4U 1820-303, although the maximum quality 
factor reached is significantly lower than for the previous 
sources. For 4U 0614+09, whereas the rise of Q is clearly 
suggested, the drop cannot be inferred. In all sources, but 
possibly 4U 1728-34, as shown at the bottom of the six 
panels, the quality factor of the upper kHz QPOs shows a 
positive correlation with frequency all the way to the highest 
detectable frequency. 

As an indication, in those files in which the QPO 
tracked is detected in at least five intervals, the mean track- 
ing timescales is ~ 230, 70, 180, 220, 140 and 290 seconds 
for 4U 1636-536, 4U 1608-522, 4U 1735-44, 4U 1728-34, 
4U 1820-303 and 4U 0614+09, respectively. These aver- 
ages are only representative because they include files in 
which the lower QPO was followed and files in which the 
upper QPO was followed by the tracking procedure, and it 
is known that the tracking timescales is longer for the up- 
per QPO than for the lower QPO, and that for the lower 
QPO the tracking timescales depend on Q; the higher Q 
the smaller the tracking timescales (see Fig 1 in Barret et 
al. 2005a). This explains why the mean tracking timescale 
for 4U0614+091 is the longest because first, there are many 
files in which the upper QPO was followed by the tracking 
procedure and second only the low Q domain of the lower 
QPO has been sampled so far with RXTE. One can also es- 
timate the minimum tracking timescale for the lower QPO. 
It is 24, 24, 48, 56, 40, and 168 seconds for 4U 1636-536, 
4U 1608-522, 4U 1735-44, 4U 1728-34, 4U 1820-303 and 
4U 0614+09 respectively. 

Using Figure 1, it is easy to determine which of the 
QPOs was followed in the data file by the tracking proce- 
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Figure 2. Quality factor versus frequency of the lower (filled 
circles) and upper QPOs (filled stars). The lower QPO has been 
sampled with a frequency bin of 25 Hz, whereas the bin is 50 Hz 
for the upper QPO. 



dure (even when a single QPO was detected). In order to 
get a better determination of the dependency of the quality 
factor of the two QPOs with frequency, we identify those 
data files in which the QPO tracked by the procedure de- 
scribed above is the lower and upper QPO respectively. The 
frequency range covered by the lower and upper QPO are 
then divided in intervals of 25 Hz and 50 Hz respectively. 
All the 8 second PDS with a tracked QPO falling into a 
frequency bin are then shifted to the central bin frequency. 
The QPOs in the resulting PDS are then fitted. The results 
of the fits are presented in Figure 2 and 3 for the quality 
factor and amplitudes of both QPOs. The trends seen in 
Figure 1 are then much clearer (this is obvious in the case of 
4U 1735-44). For 4U 1820-303, only the decaying part of the 
quality factor versus frequency curve has been sampled. The 
behavior of the upper QPO is strikingly different from the 
lower QPO one. In all, but 4U 1728-34, the quality factor 
increases all the way to the highest detectable frequencies. 
In 4U 1728-34, the peak of the quality factor is reached at 
~ 1050 Hz, beyond which a significant decrease is observed. 
In all sources, for the lower QPO, there is a clear trend for 
the amplitude to increase, then to decrease with frequency. 
The tendency for the amplitude of the upper QPO to de- 
crease with frequency is also present, indicating that the 
quality factor and amplitude of the upper QPO anticorre- 
lates towards its high-frequency end. Comparing Figure 2 
and 3, a positive correlation between the quality factor and 
amplitude of the lower QPO can be inferred during the rise 
of its quality factor. 

In Figure 4, we show the frequency difference between 
the two QPOs against the lower QPO frequency. The fre- 
quency at which the measured quality factor reached its 
maximum is also indicated, together with the spin frequency 
or half the spin frequency when available. This figure indi- 
cates that the peak of the quality factor of the lower QPO is 
not connected to the frequency difference between the two 
QPOs being equal to the spin frequency of the neutron star. 
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Figure 4. Frequency difference versus lower QPO frequency. 
When available the inferred spin frequency or half the spin fre- 
quency is shown with a dashed line. The frequency at which the 
quality factor of the lower QPO frequency reaches its maximum 
value is also shown as a vertical line. 



3 DISCUSSION 

All current QPO models face challenges in explaining the 
high quality factor observed, its dependence on frequency, 
and the significantly different behaviours of the quality fac- 
tors of the upper and lower QPOs as reported here. It is 
therefore difficult to be certain of the implications of our 
results. However, as we discussed in Barret et al. (2005a,b), 
the evidence in 4U 1636-536 for a frequency ceiling inde- 
pendent of count rate and a sharp drop in the quality factor 
at a fixed frequency is consistent with prior expectations for 
effects due to the ISCO. Such an inference would be impor- 
tant enough that it is essential to determine the implications 
in terms of neutron star masses, to examine alternate inter- 



pretations, and also to determine if the current data are at 
least somewhat quantitatively consistent with simple models 
involving the ISCO. 

3.1 Implication for neutron star masses 

Under the ISCO interpretation one would expect different 
sources to show different frequencies of peak Q and different 
frequencies at which Q extrapolates to zero (i-iower.Q^o)- Us- 
ing Slower, q=o to infer the ISCO frequencies, on should then 
verify that the latter are compatible with our understanding 
of neutron star masses and evolution. 

Figures 2 shows that the first two criteria are indeed 
met. Where it can be determined, the peak Q frequency 
ranges from ~ 750 Hz in 4U 1735-44 to ~ 900 Hz in 
4U 1728-34. Similarly, ^i OW or,Q=o ranges from ~ 850 Hz in 
4U 1735-44 to - 950 Hz in 4U 1728-34. 

Inferring the ISCO orbital frequency requires a some- 
what more specific model. All models of neutron star kHz 
QPOs that can accommodate the data suggest that the 
upper kHz QPO is close to an orbital frequency (e.g., 
Miller et al. 1998; Lamb & Miller 2001, 2003) or a vertical 
epicyclic frequency (e.g., Stella & Vietri 1998; Stella, Vietri, 
& Morsink 1999; Psaltis & Norman 2000; Abramowicz et al. 
2003; Lee, Abramowicz, & Kluzniak 2004; Bursa et al. 2004; 
Kluzniak & Abramowicz 2005), which deviates from an or- 
bital frequency by at most a few Hertz around a neutron 
star (Markovic 2000). In such models, we can use the ob- 
servation that the separation between the upper and lower 
peaks is close to either the spin frequency or half the spin 
frequency (e.g., see the panels for 4U 1636-53, 4U 1608-52, 
and 4U 1728-34 in Figure 4) to estimate the orbital fre- 
quency at the ISCO. If we simply estimate visco by adding 
half the spin frequency or the spin frequency to t'lowcr.Q^o, 
as appropriate, we find v\sco ~ 1220 Hz for 4U 1636-53, 
Hsco ~ 1230 Hz for 4U 1608-52, and ^ SC o « 1310 Hz 
for 4U 1728-34, the three of our six sources with a known 
spin frequency. From these we can estimate the gravitational 
masses of the neutron stars (see Miller et al. 1998): 



M w 2.2M (1OOO Hz/Visco)(l + 0.75j) 



(1) 



where j = cJ/(GM 2 ) ~ 0.1 — 0.2 is the dimensionless an- 
gular momentum of the star. The inferred masses of these 
sources are therefore in the range of ~ 1.8 — 2.1 Mq, which 
is higher than inferred from radio timing of double neutron 
star binaries (see, e.g., Cordes et al. 2004 for a recent review) 
but is consistent with the higher masses derived for Vela X- 
1 (Quaintrell et al. 2003) and some radio pulsars in bina- 
ries with detached low-mass stars (Nice et al. 2005). These 
masses are also compatible with realistic modern equations 
of state, which predict maximum masses for slowly rotating 
(j <C 1) stars of ~ 1.8 — 2.3 M (Akmal, Pandharipande, 
& Ravenhall 1998; Lattimer & Prakash 2001; Klahn et al. 
2006). To this level, then, identification of the rapid drop in 
quality factor as being caused by an approach to the ISCO 
is consistent with observations and theoretical knowledge 
about neutron star structure. 



3.2 Role of the magnetic field 

There do not currently exist any specific suggestions about 
mechanisms not involving the ISCO that would produce the 
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Figure 5. (top panel) The frequency ceiling of the lower kHz 
QPO of 4U 1636-536. (bottom panel) The behaviour of the qual- 
ity factor in two count rate regimes, delimited by the mean count 
rate over the whole observing set. The count rates are normal- 
ized per RXTE/PCA/PCU units and arc estimated over the full 
energy range covered by the PCA. 



observed drop in Q. Nonetheless, we can follow Barret et 
al. (2005a,b) in speculating that interaction with the stellar 
magnetic field may have an effect on the QPO coherence. 
This speculation may be supported by the fact that the 
magnetic field is likely to be involved in the QPO gener- 
ation, through the role played by the spin in the setting of 
the QPO frequencies (e.g in the forced resonance model of 
Lee et al. (2004), the magnetic field may provide the nec- 
essary periodic perturbation of the disk at the neutron star 
spin frequency). The potentially complex magnetic field ge- 
ometry, especially near the star, means that it is difficult to 
establish firm predictions about the behaviour of the quality 
factor. However, one would expect that, because quantities 
such as the Alfven radius and thickness of the disk depend 
on mass accretion rate, the frequency at which the quality 
factor is maximal and the frequency at which it extrapo- 
lates to zero would also depend on the mass accretion rate 
if the magnetosphere is involved. Evaluation of this expec- 
tation against the data is difficult because there is no direct 
measure of the mass accretion rate. 

However, it is already worth noting that the quality 
factor is not well correlated with the X-ray count rate, which 
may be considered as a proxy for the mass accretion rate, 
as shown in Figure 5 in the case of 4U 1636-536 (the same 
behaviour is observed for other sources, Barret et al. 2005b). 
Although this does not absolutely rule out interpretations 
based on plasma physics, the fundamental dependence on 
the frequency is more naturally understood if it is based on 
spacetime properties such as the ISCO. 



3.3 Advection based toy model 

In this section we describe a toy model that suggests how 
the quality factor changes as the region of the QPO genera- 
tion approaches the ISCO. Generically, we suppose that the 



frequency of the QPOs is determined by some mechanism 
in an active region (at r or b) of finite extent (width Ar or b) 
in the disk. Note that on energetic grounds the emission 
itself is likely to be produced at the stellar surface, even 
if the frequency is determined in the disk. The oscillating 
phenomenon has a characteristic lifetime (i.e. lasts a certain 
number of cycles, N cyc i es ). It is further assumed to be linked 
to the accreting matter, hence it is advected through the ac- 
tive region with a radial velocity v r . As we discuss in more 
detail below, there are then three basic parts of the observed 
width of the QPO: (1) a contribution from the finite extent 
Ar or b, (2) a contribution from the radial drift Ar drift during 
the lifetime of the oscillation, and (3) a contribution from 
the finite lifetime itself. 

We now present a simplified model for these contribu- 
tions. It is worth stressing that the model discussed below 
does not provide a physical explanation for the quality factor 
of QPOs, it simply shows how, under some reasonable as- 
sumptions, the above contributions may combine to give its 
observed frequency dependence, as found in several sources. 
Our prime focus in this paper is the sharp drop of Q for the 
lower peak at its highest observed frequencies. We there- 
fore concentrate on this first, showing that the flattening of 
specific angular momentum curve close to but outside the 
ISCO naturally leads to the observed high-frequency be- 
haviour, even in our simplified mathematical formulation. 
We then show that even if the contributions to the QPO 
widths from drift and finite lifetime are kept fixed, simple 
and physically reasonable modifications to just the contri- 
bution from A7- orb are able to reproduce both the quality 
factor behaviour of the upper peak frequency, and the lower- 
frequency behaviour of the lower peak. Therefore, although 
our model is only illustrative, it demonstrates that (1) ap- 
proach to the ISCO inevitably leads to a sharp drop in Q 
in a wide variety of models, and (2) the general picture of 
QPOs proposed by multiple authors is also consistent with 
the full frequency dependence of the quality factor of both 
peaks. 

3. 3. 1 Modeling the drop of the quality factor of the lower 
QPO 

Suppose that the lower peak is produced by the inter- 
action of either an orbital (Miller et al. 1998; Lamb & 
Miller 2001, 2003, at ^ or b) or a vertical epicyclic (Stella 
& Vietri 1998; Stella et al. 1999; Psaltis & Norman 2000; 
Abramowicz et al. 2003; Lee et al. 2004; Bursa et al. 2004; 
Kluzniak & Abramowicz 2005, at u VCTt ~ ^ or b) frequency 
at r or b with the neutron star spin. In this picture, Slower ~ 

(u olb Or ^vort) - (^spin Or ^ s pi n /2) W V OTh - (^ spin Or !/ spin /2). 

If we assume that the interaction with the spin is res- 
onant, the resonance is not expected to add significantly to 
the breadth, thus the frequency width of the lower peak is 
influenced by the following factors: 

(i) Ardrift accounting for the change of frequency over the 
lifetime of the oscillating phenomenon as the gas advects in 

(ii) A^ or b corresponding to the range of available orbital 
frequencies within the active region 

(iii) Ai^Hfe produced by the finite lifetime of the oscillating 
phenomenon 

As the true distribution of frequencies within A^ or b and 
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the true waveform of the oscillations are unknown (leading 
to uncertainties also on Afdrift), we approximate the total 
frequency width with the quadratic sum of the three contri- 
butions: 



A^total = v/(A^ dl . ift ) 2 + (A^ orb ) 2 + (A^ifc) 2 



(2) 



and the quality factor is Qi owcr = ^i ower /A!/ t otai- In the fol- 
lowing, we further assume that there are no systematics in 
our effort to track the QPO, so that the shift-and-add tech- 
nique retrieves the intrinsic width Aftotal of the QPO. We 
now estimate each of these contributions in turn. 

Estimate of Adrift — In order to produce a quantitative 
model without too many free parameters, let us consider test 
particle orbits in a Schwarzschild spacetime; we expect that 
spacetimes with moderate spin (j < 0.3) will have similar 
behaviour. In such a spacetime, the specific angular momen- 
tum of a fluid element is £ = Af 1//2 r(r — 3M) _1//2 , where we 
use geometrised units in which G = c = 1. Loss of angu- 
lar momentum will therefore lead to inward radial motion 
and thus a change in the orbit. In particular, if the orbital 
period P or b = 2n(r 3 /M) 1 ^ 2 is roughly constant during the 
lifetime of the oscillating phenomenon (as it must be for 
the observed Q > 10), then for an average specific angular 
momentum loss rate (£) the total radial drift during N cyc \ e 
cycles is 

Ar drift = 2ivN cyclc (r [i /M) 1 / 2 r 

= 2^ cyclc( rVA/) 1/2 f0fS)W. 



(3) 



As expected, close to the ISCO (r — 6M for a Schwarzschild 
geometry) where the specific angular momentum decreases 
slowly with decreasing radius, loss of angular momentum 
produces rapid radial drift. The contribution of this drift to 
the width of the QPO is then 



3 Ardrift 

Ardrift ~ — Vorb 

2 r or b 



(4) 



Estimate of A^ or b — The width A^ or b is set by the radial 
extent of the active region; this is the most difficult param- 
eter to estimate in this toy model. We will first rely on a 
somewhat specific model to get a first estimate of A^ or b, 
then we will discuss the implications of our assumption. 

For the model specific estimate we will assume that 
the active region is located at the inner edge of the disk 
where the radial velocity of the accreting gas increases dra- 
matically, for example in response to an efficient removal of 
angular momentum by the stellar radiation field or stellar 
magnetic field (e.g. Miller et al. 1998). We will assume that 
the extent of the active region corresponds to an increase 
of the radial speed comparable to the initial radial speed. 
With this assumption, we obtain: 



Ar orb 



v r (r rb) (r~6M)M 1/2 



dl/dr 2(r~3M)3/ 2 ' 
and the corresponding frequency width is 
3 Ar orb 

A^ or b ~ — ^orb 

^ r or b 



(5) 



(6) 



Estimate of Aum e — The final contribution is simply 
A^Hfe = Slower /(Tricycle). In the following we will assume 
that the number of cycles is constant over the frequency 
range considered. 
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Figure 6. Measured quality factor of the lower and upper kHz 
QPOs of 4U 1636-536, compared to the quality factor estimated 
with the toy model for the lower (solid line, see section ^. 3. II and 
for the upper (dashed line, see section [3.3.2l QPOs. 



Comparison with data — This model of quality fac- 
tor versus frequency depends on four parameters: N cyc \ e , 
v r (r xb)/(dl/dr), (£), and fisco- Assuming for simplicity 
that all these parameters are constant for a given source, 
we show in Figure 6 a comparison of a by-eye fit to the data 
for 4U 1636-536 with parameter values of 100, 0.8, 3 x 10" 8 , 
and 1220 Hz. Note that the value of {£), although apparently 
small, simply corresponds to the value that would produce 
a ~1% radial drift during the 100 cycle lifetime of the QPO. 

For these reasonable parameter values, the toy model 
reproduces quantitatively the abrupt drop of the quality fac- 
tor at high frequencies. This is the region where the main 
contribution to the QPO width is Ardrift- The simple toy 
model also reproduces the peak in the quality factor. These 
are generic features of the flattening of specific angular mo- 
mentum outside the ISCO: the active region narrows, and 
the drift speed increases, closer to the ISCO. The observed 
behaviour at high frequencies is therefore well-matched by 
these generic expectations. 

At lower frequencies, however, Q for the lower peak is 
much lower than expected in this highly simplified model. 
In addition, the quality factor for the upper peak is typ- 
ically much lower than for the lower peak when both are 
seen together. Although these are not the main focus of our 
analysis, we now evaluate whether plausible extensions to 
our picture can also match these other behaviours quantita- 
tively. 



3.3.2 Modification of the toy model 

As a guide to the necessary modifications, we show in Fig- 
ure 7 the relative contribution of the different factors to 
the total width of the QPO. This figure shows that in our 
model, A^ rb is the dominant contributor to the width for 
Slower < 800 Hz. In order to match the data, Ai/ or b should 
keep increasing with decreasing frequency (it should be ~ 60 
at 600 Hz whereas our model predicts ~ 7). 
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Figure 7. Different contributions to the total QPO width (solid 
line) estimated in the toy model, compared with the measured 
width for the lower QPO of 4U 1636-536. The three contributions 
are shown: Av^iia (long dashed line), accounting for the change 
of frequency over the lifetime of the oscillating phenomenon as the 
gas advects inwards; A^ orb (short dashed line), corresponding to 
the range of available orbital frequencies within the active region; 
and Ai^ife (dotted line), produced by the finite lifetime of the 
oscillating phenomenon. 

Therefore, in the framework of the toy model discussed 
above, we can keep the same parameters Adrift and A^nf c , 
but adjust Ai/ or b, so that a good match to the data is ob- 
tained. This is shown in Figure 8 where the difference be- 
tween the observed data points and the contribution of Ai^jf e 
and Az/drift (i-e. Ai/ or b) to the width has been approximated 
by a power law. From Af or b, one can derive an estimate from 
the data of the extent of the active region (i.e. Ar or b/r rb)- 
The bottom of Figure 8 indicates that the size of the ac- 
tive region, which determines the shape of the Q-frequency 
curve at the low frequency end, is relatively broad at low 
frequencies but decreases rapidly at higher frequencies. 

With the model described above, it is tempting to in- 
vestigate which modification would be needed to account for 
the very different behavior of the upper peak. For the upper 
peak we note that Q tends to be much lower than for the 
lower peak, that for 4U 1636-536 there is an abrupt rise in 
Qupper at the highest frequencies, and that there is no visi- 
ble drop. Remarkably, all these features can be reproduced 
just by modifying our simple model by 

Al^orb, upper = 3-upper Aforbjlowcr (7) 

where X upper IS 9j free parameter that is assumed indepen- 
dent of radius for simplicity, and the other contributions 
are kept fixed. As Figure 6 shows, :r upper = 20 matches the 
data for 4U 1636-536 surprisingly well, where Supper is com- 
puted for a given ^i OW or using the observed frequency separa- 
tion. Note that the increased breadth means that Adrift and 
Ai>iife only become important closer to the ISCO, leading to 
a sharper drop in Q that would be difficult to detect due to 
the low amplitude of the oscillation and narrow frequency 
range in which it is visible. 

Is it plausible that Af or b would behave as needed to 



Figure 8. Same as Figure 7, but with Ai/ or b fitted by a power law 
to account for the rise of the QPO width towards lower frequency 
(see section 3.3.2) (top). Ar orb /r orb as derived from Au olh fitted 
above (bottom). 

produce the observed quality factors? Although no theory 
currently exists to provide quantitative predictions, we can 
argue that this behaviour is reasonable. For example, the 
upper peak is known to have different properties than the 
lower (e.g., see Barret et al. 2005a), so additional broaden- 
ing mechanisms are to be expected. One specific proposal, 
made by Lamb & Miller (2001), is that if the upper peak 
is produced by the movement of the footprints of gas ac- 
cretion onto the neutron star, then a contribution to the 
breadth and centroid frequency of the QPO is the total az- 
imuthal angle (and its time derivative) traveled by gas from 
the active region to the stellar surface. This could be a major 
source of broadening, which would plausibly dominate the 
total width and would be proportional to Ar orb . Similarly, 
if the width of the active region is determined by angular 
momentum removal by radiation or magnetic fields, then 
regions more distant from the star could easily be screened 
more effectively, leading to a broader region and thus lower 
quality factor QPOs. Therefore, although we must wait for 
full numerical simulations to get first-principles information 
about neutron star QPO quality factors, the qualitative be- 
haviour of both QPOs is understandable in the current pic- 
tures as long as the high-frequency drop in Q is related to 
the ISCO. 



4 CONCLUSIONS 

We conclude that current RXTE data are qualitatively and 
quantitatively consistent with an ISCO-induced sharp drop 
in quality factor in several low luminosity neutron star bina- 
ries. If this is the correct interpretation, this is a signature of 
an effect unique to strong gravity, and it also provides strong 
constraints on models of matter beyond nuclear density by 
requiring that neutron star gravitational masses can exceed 
2 Mq . The most important observational work now remain- 
ing is to determine whether models not involving the ISCO 
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are plausible. The next step is thus to study the quality fac- 
tor against estimates of the mass rate such as count rate or 
hard or soft color, and more generally spectral parameters. 
This will be the scope of a forthcoming paper, using the 
unique properties of 4U 1636-536, which stands currently 
as the best candidate for the ISCO signature. 
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